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We report measurements and theory for long-range intensity correlations near the Anderson tran- 
sition of classical waves in a three-dimensional disordered material. Our ultrasonic experiments are 
designed to unambiguously detect a recently predicted infinite-range Co contribution, due to local 
density of states fluctuations near the source. We find that these Co correlations, in addition to 
C2 and C3 contributions, are significantly enhanced near mobility edges. Separate measurements of 
the inverse participation ratio reveal a link between Co and the anomalous dimension A2, implying 
that Co may also be used to explore the critical regime of the Anderson transition. 
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Wave transport in strongly scattering disordered sys- 
tems continues to reveal remarkable phenomena that 
challenge our understanding of wave physics. In address- 
ing these challenges, much can be learned by investigat- 
ing the statistics and correlations of the randomly fluctu- 
ating transmitted intensity, or speckle pattern. Intensity 
correlations of different kinds have been predicted and 
observed for classical waves under various conditions: ei- 
ther short or long range (Ci, C2, or C3), and in either 
frequency, position or angle [IHZ]- Also, correlations in 
polarization [8] or phase [9] have been investigated theo- 
retically and experimentally. For some time, there have 
been speculations about a new kind of "Co" correlation in 
the speckle pattern, with infinite range, caused by scat- 
tering in the proximity of the source [TO] [TT]. It was 
later shown that this speckle correlation is closely related 
to fluctuations in the local density of states (LDOS) of 
the waves in a random medium [T2] . This quantity is 
very actively studied both theoretically [13] and experi- 
mentally [T41IT7] . Around the mobility edge (ME) of a 
three-dimensional (3D) disordered sample, the normal- 
ized variance of LDOS fluctuations is predicted to be of 
order unity. They are also expected to become dependent 
on the sample size at, and below, the ME, due to the 
presence of states with localized or fractal spatial struc- 
ture. Recent theoretical studies even report the variance 
of LDOS fluctuations to behave critically, i.e. as a one- 
parameter scaling function of sample size and localization 
length [18 , 19 , which means they constitute a new tool to 
provide insight into the Anderson transition. However, 
despite the profusion of results concerning LDOS fluctu- 
ations, no direct experimental evidence of Co intensity 
correlations has been reported so far [20] . 

In this Letter, we present the first direct experimental 
evidence of infinite-range (Co) spatial and frequency cor- 
relations of the intensity near the Anderson transition 
of a disordered, strongly scattering 3D material. The 
experiments were performed using ultrasonic techniques 
on samples in which 3D Anderson localization of ultra- 



sound has been demonstrated previously [21] . Compari- 
son of experiment with theory, coupled with complemen- 
tary measurements designed to suppress infinite-range 
correlations when desired, allows the Co contribution to 
the correlations to be clearly separated from the other 
contributions (Ci, C2, or C3), unambiguously revealing 
the presence of large infinite-range correlations. We ob- 
serve that these correlations grow dramatically near the 
expected mobility edges in our samples. Motivated by 
the prediction that the LDOS fluctuations are closely 
related to multifractality of the wave functions through 
the q = 2 generalized inverse participation ratio (gIPR) 
[T5] H9] . we measure the anomalous dimension A2 for 
our samples in independent experiments and find good 
correspondence between this quantity and measured Co 
correlations. This clearly demonstrates the link between 
multifractality, Co, and the LDOS fluctuations. 

The samples investigated are disordered elastic net- 
works of aluminum beads, weakly brazed together to 
form slabs [21]. The monodisperse beads are 4 mm in 
diameter, and the samples have 45% porosity, consistent 
with random loose packing. This porous mesoscale struc- 
ture leads to very strong scattering with low absorption 
in the frequency range investigated (~ 0.5-2.5 MHz), a 
crucial feature for the observation of 3D Anderson lo- 
calization of ultrasound in this material. The mesoscale 
structure also leads to high contrast in the density of 
states of the aluminum matrix compared to that of the 
pores — yet another reason for anticipating strong fluctu- 
ations of the LDOS. The samples were waterproofed so 
that the experiments could be performed in a water tank 
with either vacuum or air in the pores, thereby ensuring 
that the detected transmitted waves had traveled only 
through the aluminum bead network. 

In our experiments, a broadband ultrasonic pulse is 
incident on the sample, and the transmitted waves are 
measured in the near field. The pulse is generated by 
a piezoelectric transducer, and is tightly focused to a 
spot approximately a wavelength in size on the input 
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FIG. 1. (color online) Spatial intensity correlations for the 
two types of experiments at 2.4 MHz. The scanned-source 
data show convincing evidence of infinite-range (Co) correla- 
tions, which are suppressed when only a single source point 
is used. Lines are theoretical fits using the values of parame- 
ters given in Table [I] The inset shows the single-source data 
on a log-log scale in order to reveal the extent to which the 
expected 1/Ar dependence is observed for C2 and C3 at in- 
termediate length scales. Data have been averaged over a 
bandwidth of 750 kHz, and the error bars are the standard 
deviations associated with the data's statistical fluctuations, 
which are observed to be inherently large near the Anderson 
transition. (The open symbols for the scanned-source data 
represent positions where the measurements are not as reli- 
able because of a smaller signal-to-noise ratio.) 



face of the sample. The transmitted pressure is detected 
by a sub- wavelength hydrophone, typically placed less 
than 1 mm from the sample surface (within one wave- 
length in water). To capture contributions to Co due to 
LDOS fluctuations at both the focal point of the incident 
wave (source point) and the detector, we scan both the 
source and detector over the surface of the sample. The 
recorded pressure fields p(r, t) are Fourier transformed to 
obtain the intensity I(t,uj) oc |p(r,cj)| 2 as a function of 
frequency for each pair of source and detector positions. 
The correlation function of intensity is calculated as 



(SI(r, to - ±{l)6I{r + Ar,co+ 
(I(r,w- ifi))(I(r + Ar,w+ifi)>' 



(1) 



where the angular brackets denote ensemble averaging 
and 51 = I — (I) is the fluctuation of intensity. En- 
semble averaging is done by scanning over many source 
and detector positions corresponding to the same Ar. 
For comparison, experiments with a single (stationary) 
point source were also performed, in which case the en- 
semble averaging was done only over all possible detector 
positions; this suppresses Co correlations due to LDOS 
fluctuations at the source. In what follows we will study 
spatial correlations C UJ (Ar) = C^(Ar, 0) and frequency 
correlations C w (fi) = (7^(0,0) separately. 

Figure [l] shows the spatial correlations measured near 
/ = uj/2tt — 2.4 MHz, the frequency at which Ander- 



son localization of elastic waves was demonstrated in this 
sample [21]. For both types of experiments, the correla- 
tions decay rapidly at small Ar due to Ci, with a slower 
decay due to C2 and C3 that extends out to Ar ~ 10 mm, 
beyond which C 0J (Ar) becomes independent of distance. 
For the data where the source position is varied, an 
asymptotic value of order unity is seen for the correla- 
tions, showing clear evidence of a Co term due to LDOS 
fluctuations at the source. By contrast, no infinite-range 
correlations are seen for the single-source data, consistent 
with the fact that the LDOS at the source position does 
not fluctuate in this case. 

To gain further insight into this behavior, we compare 
our experimental data with theoretical calculations. We 
compute Ci, C2, C3, and Co correlation functions assum- 
ing weak disorder and write the full correlation C UJ (Ar) 
as a function of three fit parameters: A, Cq^ , and Cq OU ^ . 
The parameter A quantifies the magnitude of C2 and 
C3 correlations, Cq 11 ^ characterizes the magnitude of the 
genuine Co correlation due to the LDOS fluctuations at 
the source point, and Cq OU ^ measures the amplitude of 
the short-range contribution to Co due to scattering in 
the vicinity of both detectors when the latter are close to 
each other [22]. Cq^ is the asymptotic value of C u (Ar) 
for Ar — >• 00. The solid lines in Fig. [I] show the re- 
sults of performing a joint weighted fit of these theoret- 
ical predictions to both the single- and scanned-source 
data at each frequency, thereby determining the values 
of the parameters shown in Table [IJ In these fits, we 
account for the fact that Cq U ^ contributes only to the 
scanned-source correlations and set Cq^ = for fit- 
ting the single-source data; also, since the detector ge- 
ometry is the same for both experiments, Cq OU ^ is con- 
strained to have a common value for the two curves. Note 
that for white-noise uncorrelated disorder and point-like 
source and detector in an infinite disordered medium, 
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7r/k£, where k is the wavenumber in 



the medium and i is the mean free path [10]. In our ex- 
periments, however, both the source and detector have 
finite extent (which may differ from each other), and the 
finite size of the aluminum beads inevitably results in 
some short-range structural correlations. Therefore, we 
expect in general that C^ ^ Cq OU ^ ^ ir/ki [TT] , Fig- 
ure [T] provides strong evidence that the large asymptotic 



value of C u (Ar — » 00) = C { 



(in) 



1 for the scanned-source 



experiment is due to Co correlations. 

Similar behavior, with Cq 11 ^ ~ 1, is observed over a 
broad frequency range from 1.6 to 2.8 MHz, where in- 
dependent measurements of the dynamic transverse con- 
finement of the transmitted intensity [21 J indicate that 
ultrasound is still localized, with similar values of the 
localization length £ (£ w L = 14.5 mm for this sam- 
ple in this frequency range). At lower frequencies, at 
least one ME must exist, since previous measurements on 



TABLE I. Fit parameters, with uncertainties in parenthe- 
ses. The uncertainties are given by the standard deviation 
of the parameters. For a point source and detector, the nor- 
malized variance, (7(0,0), depends on all three parameters 
as C(0, 0) = 1 + 2[A + C<j in) + C ( ° ut) }. By contrast, the in- 
finite range contributions depend independently on the dif- 
ferent contributions to Co, with the asymptotic values of 
the scanned-source C(Ar, 0), the single-source (7(0, Q) and 



the scanned-source (7(0, Q) being equal to C<j in) , C^ outj , and 



-v(out) 



Cq 1U ^ + Cq 011 ^ , respectively. 



Parameter 2.4 MHz 0.97 MHz 1.07 MHz 1.11 MHz 



Spatial correlations 



A (single) 
A (scanned) 
C, 



(in) 


(out) 



0.50(0.02) 0.48(0.02) 1.29(0.04) 1.59(0.06) 

2(1) 0.8(0.2) 6(1) 6(3) 

1.3(0.2) 0.42(0.02) 1.06(0.08) 7.8(0.5) 

0.4(0.2) 0.8(0.2) 1.5(0.4) 7(1) 



Frequency correlations 



A (single) 
A (scanned) 



^(out) 
°0 



Q T h/27T (kHz) 



0.2(0.2) 
0.7(0.1) 
0.32 (0.03) 
0.62 (0.06) 
4.35 (0.09) 



0.8(0.3) 
5.2(0.8) 
0.9(0.3) 
1.3(0.1) 
7.2(0.1) 



these samples revealed that ultrasound propagates diffu- 
sively at the much lower frequency of 200 kHz [21 . To 
investigate the long-range correlations as a ME is ap- 
proached, experiments were performed at intermediate 
frequencies, between these well established diffusive and 
localized regimes [21], with representative data at fre- 
quencies near 1 MHz being presented in Fig. [2] Both 
spatial and frequency correlations increase significantly 
with frequency when a ME, which we estimate to be at 
approximately 1.1 MHz, is approached. In particular, 
the asymptotic value of the scanned-source spatial cor- 
relations, Cq U \ increases from 0.4 to almost 8 over the 
range of frequencies illustrated in Fig. [2|a). 

The frequency correlations also show large increases in 
Co over this frequency range [see Fig. j2^b)] [23]. Cu(Q) 
contains infinite-range contributions from scattering both 
near the source and near the detector, i.e. both 
and CQ OUt ^ contribute to the asymptotic value of 0^(0) 
for large Q [22]. The single-source measurements (which 



suppress Cq^) show that Cq OUT,) increases from 0.6 to 1.3 
between 0.97 and 1.07 MHz. By comparing the best-fit 
values (see Table [i]), we see that for the scanned-source 
case, C<j in) and Q out) are of the same order of magni- 
tude, as could be expected from the roughly symmetric 
arrangement of the experiment [24] . 

The C2 and C3 correlations, quantified by the param- 
eter A, also increase with frequency around 1 MHz, as 
found from the data for both spatial and frequency cor- 
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FIG. 2. (color online) Spatial (a) and frequency (b) corre- 
lations measured near 1 MHz, showing the increase of long- 
range correlations near a mobility edge. The inset shows the 
characteristic 1/v^ behavior expected for C2 and C3 corre- 
lations. For both plots, data have been averaged over a band- 
width of 25 kHz, except at the highest frequency (1.11 MHz), 
where the data are changing too rapidly with frequency to be 
meaningfully averaged. These rapid variations with frequency 
near 1.11 MHz also complicate measurements of frequency 
correlations, which are therefore not shown here. The error 
bars are calculated as in Fig. ^ Lines show the fits using the 
parameters given in Table [I] 



relations (see Fig. [2] and Table]]]). Because A oc l/(k£*) 2 
to leading order [22] , where £* is the transport mean free 
path, the increase of A corresponds to a decrease of k£* as 
the ME is approached. In addition, the values of A found 
from the fits are always larger in the scanned-source case. 
Keeping the source fixed not only suppresses c£ n) but 
also reduces the magnitude of long-range C2 and C3 cor- 
relations because the latter contain contributions from 
scattering in the vicinity of the source. This effect does 
not preclude the clear identification of Co that stands out 
by its infinite range in both space and frequency. 

The frequency dependence of the asymptotic value of 
the spatial intensity correlation function between 0.6 and 
1.4 MHz is shown in Fig. |3ja). These data are the av- 
erage of the measured correlations for Ar between 25 
and 50 mm, where C L0 (Ar) is found to be indepen- 
dent of distance, providing accurate measurements of 
Co(Ar — » 00) = C^ when the source is scanned. It in- 



4 



+ 



< 



10 



1 

-0.2 
-0.4 
-0.6 
-0.8 





o 



1 10 " 

§ 10 _: 
?— i 

H 



► Scanned source 
" ° Single source 


' 1 ' 1 


1 1 1 1 1 1 


i ' ^ 






■ 1 i 1 i 1 i 1 i 

' vs ~n 

V 


(b)t 

\ ^ ) 


1 1 1 1 1 1 


1 1 1 


— i — 1 — i — 1 — i — F 

~ (en 


ft Vi.tT. 


v/t . .Tli 



0.6 0.8 1.0 1.2 

Frequency (MHz) 



1.4 



FIG. 3. (color online) (a) Frequency dependence of the 
asymptotic value of the spatial correlation Co;(Ar — > oo) = 
Cq^ for single and scanned point sources, (b) anomalous di- 
mension A2 of the gIPR for q = 2, and (c) the amplitude 
transmission coefficient. The arrows in (c) indicate the reso- 
nant frequencies of individual aluminum beads. Cq 111 ^ and A 2 
show large increases in magnitude near the upper band edges, 
where we expect mobility edges to occur. The grey vertical 
line indicates the location of the ME near 1.1 MHz, estimated 
from transverse confinement measurements . 



creases rapidly with frequency near 0.78 and 1.11 MHz, 
reaching values up to 13 here, and even as high as 30 in 
other experiments — by far the largest values of Co ever 
reported. Comparison of these results with the ampli- 
tude transmission coefficient [Fig. |3jc)] reveals that the 
frequencies where Co increases rapidly coincide with the 
upper edges of pass bands in these disordered structures. 
In the band gaps, the transmission becomes too small to 
yield measurable long range correlations. As explained in 
Refs. [21j[25], these band gaps are not due to Bragg scat- 
tering, as in phononic crystals [26]. Instead, they arise 
between pass bands formed from coupled resonances of 
the beads when the coupling is sufficiently weak. 

Near the upper edges of the pass bands, where the av- 
erage density of states decreases, mobility edges between 
extended and localized states may be expected [27]. Ev- 
idence that mobility edges do indeed occur near these 
band edges has been obtained through separate measure- 
ments of increased spatial confinement of the transmit- 
ted intensity near the upper band edges relative to the 
pass band centers, using the method developed by Hu 
et al [21]. This evidence is most compelling for the 
ME near 1.1 MHz, which is indicated by the vertical line 
in Fig. [3] Additional evidence can be inferred from the 
large increases that are found in the normalized intensity 
variance, C 0J (Ar = 0, Q = 0), near the upper band edges 



(e.g. Fig. [2] and Ref. [21 J). Thus, the large increases in 
Co near 0.78 and 1.1 MHz must be due to large LDOS 
fluctuations near Anderson transitions in these samples, 
suggesting that Co is sensitive to critical effects. 

This interpretation of the striking increase in Co near 
the band edges is further supported by measurements of 
the anomalous multifractal dimension A2, which char- 
acterizes the length-scale dependence of the inverse par- 
ticipation ratio (IPR) P 2 ~ L~ d ~ A2 [28 . The sig- 
nificant decrease in A 2 with frequency near the upper 
band edges [Fig. [3^b)] is consistent with the predicted 
behavior near the Anderson transition, where A 2 should 
become increasingly negative, varying from in the dif- 
fuse regime to —2 deep in localized regime [29]. Since the 
source and detector in our experiments are point-like, it 
is likely that a single mode dominates at any frequency, 
so we expect the IPR calculated from the intensity J(r) 
and from the LDOS p(r) to be equal p~5j [19]. Then, 
P 2 = L- d (p 2 )/(p) 2 = L- d [C (oo) + 1], and we predict 
that log [Co (00) + 1] oc — A 2 . Within experimental er- 
ror, the frequency dependencies of C^Ar — >• 00) = C^ 
and A 2 [Figs. j3^a) and (b)] are consistent with this pre- 
diction. Thus, not only do the infinite-range correlations 
and the IPR show evidence of transitions from extended 
to localized behavior near the upper band edges, but the 
correspondence between these measurements verifies the 
link between Co, LDOS fluctuations, and multifractality 
experimentally. 

In conclusion, infinite-range intensity correlations have 
been measured directly in a strongly scattering 3D 
"mesoglass" for which Anderson localization of ultra- 
sound was unambiguously demonstrated previously [21] . 
Measurements are consistent with the diagrammatic the- 
ory when large magnitudes of both long-range (C 2 and 
C3) and infinite-range (Co) terms are assumed. By vary- 
ing the ultrasonic frequency, we have been able to in- 
vestigate the growth not only of C 2 and C3 but also of 
Co correlations as the Anderson transition is approached. 
Infinite-range correlations of order unity are found over 
a broad range of frequencies, reflecting the high LDOS 
contrast that can be achieved in our samples. The mag- 
nitude of these Co correlations is seen to increase dramat- 
ically near a ME. These Co results are mirrored by the 
frequency dependence of the anomalous dimension A 2 , 
which characterizes the spatial structure of the inverse 
participation ratio. Our independent measurements of 
these two quantities establish a link between Co and A 2 , 
revealing that Co can be used to probe the Anderson 
transition. The possibility of exploiting our findings to 
experimentally investigate critical behavior at the Ander- 
son transition, by focusing on the possible one-parameter 
scaling of Co near the ME, is a promising new avenue for 
future research. 

This work was supported by NSERC and by a PICS 
program of CNRS. 
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Supplemental material 



Long-range correlation C2 



INTRODUCTION 

This document provides the summary of theoretical 
results that we used to fit the data in the main text of 
the paper. We compute the short-, long- and infinite- 
range correlation functions of intensity under assumption 
of weak disorder k£ ^> 1 , where k is the wave number and 
£ is the mean free path due to disorder. 



SPATIAL CORRELATIONS 
Definitions 

We consider the spatial correlation function of inten- 
sity fluctuations 5J(r, a;) = J(r, a;) — (J(r, u;)): 



CUr,r') 



(£I(r,^I(r») 
(/(r,a,))</(r',a,)>- 



(SI) 



To lighten the notation, we will omit the subscript 'a;' 
from here on, keeping in mind that all measurements 
are performed for waves at the same frequency uo. Typ- 
ically, C(r, r') is a decaying function of Ar = |r — r'| 
with C(r, r' = r) = (SI(r) 2 ) / (I(r)) 2 being the normal- 
ized variance of intensity, and C(r, r') = in the absence 
of intensity correlations. It is convenient to split C(r, r') 
in several parts: C = C\ + C2 + C3 + Co, each of C$ 
originating from different physical processes [SlJ [S2] . 



Short-range correlation C\ 

In the bulk of a disordered medium and far from 
boundaries, the short-range contribution to C is |S1[ IS3] 



sin&Ar\ 2 

J exp(- 



&Ar 



-Ar/*), (S2) 



where fc = 27r/A, £ is the scattering mean free path. At 
the surface of a disordered sample, the spatial correlation 
C\ is modified due to the anisotropic angular distribution 
of intensity [S4 : 



C^ rface (r,rO 



1 



sin/cAr J\(kAr) 



IA + ^ 

x exp(— Ar/£) 



kAr kAr 
= h(kAr, M), 



(S3) 



where A = zq/£* with zo the extrapolation length en- 
tering the boundary conditions for the average intensity, 
and we defined the function h(kAr, k£) that will be used 
in the following. This result is largely independent of the 
spatial extent of the source (plane wave, beam of finite 
size or point source) and has been tested experimentally 



In contrast to Ci, the long-range contribution C2 de- 
pends on the spatial extent of the source. It is not easy 
to calculate for an arbitrary source. In addition, we have 
to make an assumption of weak disorder (k£ ^> 1) to 
compute the diagrams corresponding to C2. 

Transmission of a plane wave through a slab. 
We assume that a slab of thickness L ^> £ and trans- 
verse extent W ^> L is illuminated by a plane wave. The 
spatial correlation of intensity is calculated at the oppo- 
site side of the slab, as a function of transverse distance 
Ar = |r-r'| > t 



a 



plane wave 



(Ar) 
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and I* is the transport mean free path. 

Equation (S4) applies for Ar ^> £* only. It is clear 
that its divergence at Ar — » is an artifact of approx- 
imations made during its derivation. Physically, we ex- 
pect C*2(Ar = 0) c± C2(Ar = £*). A more precise shape 
of C2(Ar) at small Ar < £* can be obtained by paying 
more attention to small q and avoiding the limit q£* <C 1 



which is tacitly taken in the derivation of Eq. (S4). We 
then obtain a longer but more accurate expression for C2 : 



C plane wave 
2 



(Ar) 



where 



F 2 (x,y) 



rC 



du 



2(H*) 2 
J (ux) 



(S6) 



(uysmhu) 2 ^ ^ ^ 
x [sinh[2u(l - y)] - 2u(l - y)] 
+ sinh 2 [ii(l — y)\ [smh(2uy) — 2uy]} 

1-y, x = 0, 
~ { y/x, y < x < 1, 



(S7) 



ex e" 



x > 1. 



A comparison of Eqs. (S4) and (S6) is shown in Fig. SI 



The latter equation, in contrast to Eq. (S4), allows us to 
obtain the value of C2 for Ar 
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Point source in the infinite medium. Even if 
this geometry seems simple, the calculation of C2 ap- 
pears quite involved. If we define the center of mass 
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FIG . SI. Comparison of Eqs. ( |S4| (shown by circles) and 
\S6\ (shown by the solid red line) for the long-range correla- 
tion function of intensity fluctuations. The dashed line shows 
t/Ar. 



R = |(ri + and the difference Ar = ri — T2 coordi- 
nates, we can obtain simple results for Ar _L R: 

3 t 



^opoint 



3 (Ar) ~ ±- ( 



(S9) 



2{kl*) 2 Ar' 

with the '+' sign for Ar <C R and the ' — ' sign for Ar ^> 
R. 

Transmission of a tightly focused beam through 
a slab. This situation is realized in our experiments and 
is somewhat intermediate with respect to the two previ- 
ous cases (the sample is a slab, but the source is point- 
like). Because the results for the plane wave incident on 
a slab ( |S6| ) and the point source in the infinite medium 
(S9) coincide for £* < Ar < L, we expect that the same 
result will also hold for the tightly focused beam. We 
will therefore use: 
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(Ar) to be different 
from both Eqs. (|S6j) and (|S9j), but because its magnitude 



is already small at such large distances, it will not play 
a significant role in the fits to the experimental data. 

Short-range part of C<i- The calculation leading to 
Eqs. (S6), (S9) and (S10) also yields short-range terms 



that are rarely mentioned but exist. The full expression 
for C2 including both long- and short-range contributions 
is 
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FIG. S2. Diagrams contributing to Co correlation function. 
r is the source position. The diagram (a) is the original 
long-range one S2 ; it is independent of Ar = |r — r'|. The 
diagram (b) is short-range. It was calculated in Ref. [S8j . 
The diagrams (c) and (d) are both short-range and were not 
considered previously. A complex conjugate diagram should 
be added to each of the diagrams. 



Long-range correlation C3 

The calculation of the spatial C3 correlation function 
for a beam focused on the surface of a 3D disordered slab 
is a complicated task that we did not succeed in accom- 
plishing. However, the structure of the result may be 
anticipated from the diagrams involved in the calcula- 
tion |Slj : we expect short- and long-range terms similar 
to C2. The magnitude of C3 is expected to be of order 
1/(H*) 4 . Hence, an approximate expression for C3 may 
be written as 
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Infinite-range correlation Co 

Similarly to C2 and C3, Co correlation also contains the 
'interesting', infinite-range part and the 'trivial', short- 
range one. The full expression is found by summing the 
diagrams of Fig. [S2) 
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(S13) 
(S14) 

(S15) 
(S16) 



where Cq^ is the genuine, infinite-range correlation that 
survives at large Ar [S2]. It results from the scattering 
near the source and is related to the variance of the lo- 
cal d ensity of s tates at ro [S9]. In contrast, the terms 



(Sll) (S14) and (S16), which are proportional to Co° ut \ result 



from the scattering near the detection points ri and r2. 
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They are appreciable only at small Ar = |ri — r 2 1 . For the 



white- noise uncorrected disorder, C, 



(in) _ ^f(out) 



In our experiment, the disorder is correlated and the sym- 
metry between the 'point-like' excitation and the 'point- 
like' detection may be broken because neither is actually 
point-like and the effective sizes of the excitation and de- 
tection areas may differ, so that Cq^ ^ Cq ^ ^ n/k£. 
Moreover, these parameters are not universal and will 
depend on the microscopic structure of the disordered 



sample |S10| . We use Cq and Cq ) as free fit param- 
eters when comparing theory to the experimental data. 

The functions fb(kAr,k£) and fd(kAr,k£) are rapidly 
decaying functions of kAr: 
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The behavior of these functions is illustrated in Fig. [S3| 



Full expression for the correlation function 



Adding up Eqs. Q, ( |STT| ), ( [S12] ) and ([ST3])-([ST6 ) 
we end up with an expression that can be used to fit 
experimental data: 
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, (S19) 



where A = [3/2(H*) 2 + const/ {kt) A ]F 2 (0, t/L). A, 
Cq 1 ^ and Co° ut ^ are the unknown fit parameters. 



Accounting for the size of the detector 

In the experiment, the acoustic field is measured 
very close to the sample surface with a disk-shaped hy- 
drophone of radius b = 0.2 mm. In order to take into 
account the size of the hydrophone, we assume that the 




FIG. S3. Functions fb and fa describing the short-range part 
of Co correlation function. 



measured quantity is not the intensity J(r) at a point r 
but the intensity averaged over a disk of radius b centered 
at r: 



J(r) = 4* / Hr f )d 2 r f 
Kb 1 J b(r) 



(S20) 



where 6(r) denotes a disk of radius b centered at r. The 
correlation function of X(r) can be then obtained from 
the correlation of J(r) by a double spatial integration: 

Cx(Ar = n - r 2 ) = * f dh[ [ d% 

I 7 ™ ) J bin) Jb{v 2 ) 

x C^Ar^r;-^). (S21) 

FREQUENCY CORRELATIONS 
Definitions 

The frequency correlation function of intensity fluctu- 
ations 6I(r,(j) = i(r,co>) — (J(r, a;)) is defined as 

r (r o\- ( 5I ( r >" + &) si(r,w-^)) 
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where we assume that the average intensity is indepen- 
dent of frequency in the frequency band under consider- 
ation: (J(r,a;+ = (I(r,u- = (I(r,u)). Once 
again, we will omit the subscript 'a;' of C from here on. 
The behavior of C(r, Q) with ft is similar to the behav- 
ior of the spatial correlation function C(Ar) with Ar: it 
decays and has both short- and long-range parts. 



Short-range correlation C\ 

The short-range part of C can be easily calculated in 
transmission of a plane wave through a slab of thickness 
LEB EH: 



Ci(ft) 



L smh\at) 
£* a£* sinh aL 



(S23) 



where aL = TTy/iQ/QTh, ^Th = tt 2 D/L 2 is the Thou- 
less frequency, D is the diffusion coefficient of the wave, 
and we neglected corrections due to boundary conditions, 
assuming L ^> I* , zq . 

For a point source at the origin in the infinite medium 
we have 



Ci(R,ft) = |exp(-a#)| 2 , 



(S24) 



with similar definitions aR = 7r^/^V^Th5 ^Th = 
tt 2 D/R 2 . 

Both correlation functions (S23) and (|S24|) oscillate 
and decay roughly exponentially with y/ft/ftxh, so that 
no correlation is left for ft ^> ftxh- We will adopt Eq. 



( S23 ) in the following. 



Long-range correlation C2 



Transmission of a plane wave through a slab. 

From a calculation following Refs. [S7, Sll we found the 
following result: 
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(S26) 
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.(S27) 



Integrations in this equation can be carried out analyt- 
ically, resulting in a long expression that we do not re- 



produce here. Then the integral in Eq. (S26) can be 
calculated numerically. 



In the limit of a thick slab L > f , Eq. jS25| > yields 
a function that depends mainly on ft/ftxh as far as 
ft/ftTh ^ 1, see Fig. [S4] (top). In the limit of large 
ft/ftxh 00 we find 
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as illustrated in Fig. |S4| (bottom). 

Point source in the infinite medium. Here we 
need to introduce a spatial cut-off ~ £* to avoid the path 
crossing (Hikami box) being closer than t* to the detec- 
tor. The results then should be understood as depending 
on the precise value of this cutoff: 
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In the limit of R ^> I* that is of interest for us here, we 
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FIG. S4. Frequency correlation of intensity fluctuations C2 
in transmission of a plane wave through a disordered slab for 
small (top) and large (bottom) values of O. 
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FIG. S5. Frequency correlation of intensity fluctuations C2 
at a distance R from a point source in the infinite medium for 
small (top) and large (bottom) values of 
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(S31) 

, ft>ft Th (S32) 



The behavior of Ci at small and large $7 is illustrated in 
Fig.|S5| 

Transmission of a tightly focused beam through 
a slab. We assume that the C2 correlation function for a 
beam focused on the surface of a disordered slab is similar 
to the one for the plane wave, except for the magnitude 
of C 2 |5TT] : 
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with F 2 (aL,l*/L) defined by Eq. (S26). 



Infinite-range correlation Co 

The frequency dependence of the Co correlation func- 
tion is obtained by calculating the diagrams of Fig. S6 
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FIG. S6. Diagrams contributing to Co correlation function, 
ro is the source position; cji,2 = db |0. The diagrams (a) 
and (b) are independent of f2 = cji —002 as far as |Q| <C cji,2. 
A complex conjugate diagram should be added to each of the 
diagrams. 



We obtain: 



C { *\tl) = C { 



(in) 



✓~<(out) 
°0 1 
(in) 



c^ d \n) = c^ mi c 1 (n) 



(out) . 



(S34) 
(S35) 
(S36) 
(S37) 
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Full expression for the correlation function 



Adding up all the contributions and assuming (as in 
the case of spatial correlations) that the behavior of C3 
as a function of frequency is similar to that of C2, we fi- 
nally find the full expression for the frequency correlation 
function: 
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where A 



c, 



(in, out) 




l/(k£- 



*\2 



Note that the constants A and 



here should be the same as in Eq. (S19) for the 



spatial correlation. 
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